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ABSTRACT 

Context. The study of nebulae around Wolf-Rayet (WR) stars gives us clues about the mass-loss history of massive stars, as well as 
about the chemical enrichment of the interstellar medium (ISM). 

Aims. This work aims to search for the observational footprints of the interactions between the ISM and stellar winds in the WR 
nebula NGC 6888 in order to understand its ionization structure, chemical composition, and kinematics. 

Methods. We have collected a set of integral field spectroscopy observations across NGC 6888, obtained with PPAK in the optical 
range performing both 2D and ID analyses. Attending to the 2D analysis in the northeast part of NGC 6888, we have generated maps 
of the extinction structure and electron density. We produced statistical frequency distributions of the radial velocity and diagnostic 
diagrams. Furthermore, we performed a thorough study of integrated spectra in nine regions over the whole nebula. 
Results. The 2D study has revealed two main behaviours. We have found that the spectra of a localized region to the southwest of this 
pointing can be represented well by shock models assuming = 1000 em~ 3 , twice solar abundances, and shock velocities from 250 to 
400 km s -1 . With the ID analysis we derived electron densities ranging from <100 to 360 cm 4 . The electron temperature varies from 
~7700 K to ~ 10 200 K. A strong variation of up to a factor 10 between different regions in the nitrogen abundance has been found: 
N/H appears lower than the solar abundance in those positions observed at the edges and very enhanced in the observed inner parts. 
Oxygen appears slightly underabundant with respect to solar value, whereas the helium abundance is found to be above it. 
We propose a scenario for the evolution of NGC 6888 to explain the features observed. This scheme consists of a structure of multiple 
shells: i) an inner and broken shell with material from the interaction between the supergiant and WR shells, presenting an overabun- 
dance in N/H and a slight underabundance in O/H; ii) an outer shell very intense in [OIII]/15007A corresponding to the main sequence 
bubble broken up as a consequence of the collision between supergiant and WR shells. Nitrogen and oxygen do not appear enhanced 
here, but helium appears enriched, iii) And finally it includes an external and faint shell that surrounds the whole nebula like a thin 
skin representing the early interaction between the winds from the main sequence star with the ISM for which typical circumstellar 
abundances are derived. 
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1. Introduction 

Wolf-Rayet (WR) stars are the evolved descendants of massive 
OB stars (M > 25Mq), and they are believed to constitute the 
previous phase in their evolution before supernova explosion. 
Many of them have been shown to be surrounded by nebular 



emission in a structure of shells (ring nebula). Johnson & Hogg 
( |1965| l were the first to suggest that these nebulae were due 
to material ejected from the star sweeping up the surrounding 
interstellar medium (I SM). |Pikel'Ner & Shche glov| ( |1969| l and 
| Weaver et al. (1977]) both made models deriving analytical 
solutions about dynamic evolution of shock bubbles created 
by interaction between the ISM and the stellar wind. One of 
the most comprehensive models describing the basic optical 
morphology of bubbles in WR stars was produced by |Garcia-| 
Segura et al. ([1996). Their hydrodynamical simulations take the 



* Based on observations collected at the Centro Astronomico 
Hispano Aleman (CAHA) at Calar Alto, operated jointly by the 
Max-Planck-Institut fur Astronomie and the Instituto de Astrofisica 
de Andalucfa (CSIC). 
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evolution and mass-loss history into account over the lifetime 
of the central star from the main sequence (MS) to the WR 
phase. In summary, stellar winds from the MS star sweep up 
the ISM, forming a hot bubble surrounded by a thin shell; later, 
dense material from the red supergiant (RSG) winds partially 
fills up the bubble creating another internal shell; finally, in the 
WR stage, the star develops another mass-loss episode with a 
fast wind that reaches the RSG shell and can eventually break 
up the MS bubble, creating a wind-blown bubble (WBB). 

NGC 6888 (S105) is a proto-typ ical WBB associated w ith 
the WN6 star WR 1 36 (HP 192163) (|Johnson & Hogg|1965) . It 
was first reported by Sharpless ( 1959[ ), and its relative proximity 
and large optical angular size has made it one of the best-studied 
examples of this class of objects, which are observed in a wide 
wavelength range. The basic parameters of star and nebula are 
summarized in Table Q] 

The ellipsoidal shell appears to be geometrically prolate 
and to have a highly filamentary structure (Lozinskaya 1970 1. 
The exp ansion velocity of the shell varies in the range 55-110 
km s _1 ( Lozinskaya| 1970 Marston & Meaburn||l988[ [Esteban 
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Table 1. Physical properties of WR 136 and NGC 6888. 





Object 


Parameter 


Value 


Kelerence 


WR136 


(a,6) (J2000) 


(20:12:06.55,+38:21:17.8) 


van der Hucht (2001 1 




Spectral type 


WN6 


van derHucht|l|2UUl| 






lb 


Hamann et al. (2UUb) 




Distance (kpc ) 


1.26 


vanderHucht|(|2U0l| 




Ra (kpc ) 


1U 


Hsteban et al. (1992) 




M v (mag) 


-4.69 


Hamann et al. (2U06) 




feb-v (mag; 


0.4b 


Hamann et al. (2UUt>) 




NGC 6888 


Angular size (arcmin z ) 


12x18 


Chuetal.||1983|l 




Vexo (kms" 1 ) 


55-110 


Lozinskaya ( 1970 1 




^ionized {M^) 


b 


Wendkeretal. (197b) 




^neutral (*?0) 


40 


Marston & Meaburn ( 1988| 





|& Vilchez||1992| among others). Spectroscopy of the nebular 
shell shows that the ionized gas is enriched with nitrogen and 
helium (Wendker e t al.||197"5) and is slightly underabundant in 



oxygen (Kwitter 1981 Esteban & Vilchez 1992 1. This could 



be explained as a consequence of the mixing of interstellar and 



stellar wind material ( Kwitter|1981| l and the transport from the 



core of the star to the shell during the RSG phase ( |Esteban 



& Vilchez||1992| l. |Marston| ( |1995| > presented IRAS images 
obtaining two concentric shells associated with WR 136. With 
the inferred masses, he proposed that material of the inner shell 
comes from the mass lost by the star in the RSG phase and 
the outer she ll formed from the early O star phase. Gruendl 
j20Q0| > observed Ha and [OIII]/15007A emission from 



et al. 



NGC 6888, finding an offset between both images and suggest 
that these offsets are due to the shock from the WR bubble 
expanding into the circumstellar medium. This contrast between 
the Ho- and [ODTJ/15007A morphology is shown better in the 
images from the WFPC2 on the HST ( [Moore et al.|2000) . 

Bubbles blown by massive stars are filled by the fast stel- 
lar wind that interacts strongly with the ISM, and some fraction 
of the kinetic energy is converted into thermal energy heating 
the gas at X-ray emitting temperatures (Wea veTet al.| T977] >. 
NGC 6888 was the first WB B that has been dete cted in X-rays, 
it was made with Einstein by Bochkarev ( 1988). The Wrigge et 
[aT7|(| 1 994 ) X-ray observation did not detect a filled bubble emit- 
ting as models predicted, but a filamentary X-ray structure was 
instead found following the Ha features of the nebula. In addi- 
tion, spectra in the shell were fit to a double -temperature plasma 
model with a high-temperature component at T ~ 8.5 x 10 6 K 
and a low -temperature component at T ~ 1.3 x 10 6 K without 
detecting significant temperature variation between different ar- 

|2005|l . In the last X-ray study 
have found a bi-modal 



eas of the nebula (Wrigge et al. 



of this nebula, Zhekov & Parkff20lT 



distribution of the X-ray emitting plasma, and made models to 
explain how the soft-X rays come from evaporation of clumps. 
They also discuss whether the heating conduction is efficient in 
the hot bubble. 

NGC 6888 represents one of the most suitable candidates 
for searching for the observational footprint of the interactions 
between ISM and stellar winds. In addition, it can give us clues 
about the mass -loss history of massive stars, as well as about 
the enrichm ent of the ISM ( |Esteban et al.|1992[ |Monreal-Ibero| 
|et al.|2010| i. Despite all the studies already done on NGC 6888, 
the spatial extent of the chemical enrichment and the homogene- 
ity of the physical properties and structure are still unknown 



, since no 2D measures of their spectroscopic properties (and 
imaging, simultaneously) have been done so far. To study the 
morphology and chemodynamics in 2D, we decided to include 
NGC 6888 as a target in our programme of integral field spec- 
troscopy (IFS) to perform the spatial analysis of spectral maps 
and compare them with ID studies of selected integrated spectra. 

The paper is organized as follows. The observations and data 
reduction are described in Sect. [2] Description, analysis, and re- 
sults of the physical conditions of the 2D study are presented in 
Sect. [3] and Sect. |4]is devoted to study selected ID spectra. The 
discussion and the summary of the main conclusions are given 
in Sects. |5]and|6j respectively . 



2. Observations and data reduction 

2.1. Observations 

Observations were carried out in 2005 July at the 3.5-m tele- 
scope of the Centra Astronomico Hispano Aleman (CAHA) at 
the observatory of Calar Alto (Almera, Spain) with the Potsdam 
Multi-Aperture Spectrograph instrument (PMAS) ( |Roth et aE] 
2003] ) in PPAK mode ( |Kelz et al.|2006) . The PPAK configura- 



tion consists of 331 science fibres with a diameter of 2.7 arcsec 
covering a hexagonal field of view (FoV) of 74x65 arcsec 2 . The 
surrounding sky is sampled by 36 additional fibres distributed in 
six bundles located following a circle at about 90" from the cen- 
tre. Furthermore, there are 15 fibres for calibration purposes (see 
Fig. 5 in |Kelz et al.|2006) . 

We chose the zones of NGC 6888 to be mapped based mainly 
on the narrow band images observed by our group in the Isaac 
Newton Telescope (INT) with the Wide Field Camera (WFC). 
Figure [T] shows an RGB composited image of the nebula; Ha 
and [OIII] emissions represented by red and green colours re- 
spectively, the [Oil] emission, in blue, is fainter but noticeable 
as the pink shade in the mix with Ha. Here, we present IFS ob- 
servations of NGC 6888 in four zones of the nebula (see Fig. [TJ: 
the first zone (1) in the northeast area, was chosen because it can 
be a good candidate for presenting signs of shocked gas, and it 
is known that this area emits in X-rays. Two other pointings, 
(2) and (3), were chosen in the outer zones of NGC 6888 where 
Ha emission is very weak but an [OIII] skin is observed in Fig. 
[T] Finally, another zone (4) was selected to study the centre of 
the nebula where we can observe a striking dark region in all 
the wavelengths or a knot that could be associated to our ob- 
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Fig. 1. RGB composited image of NGC 6888 taken with the WFC at the Isaac Newton Telescope (INT). Red corresponds to the Ha 
emission, green to [OIII], and blue to [Oil]. [Oil] emission is very faint but noticeable as the pink shade in the mix of Ha and [Oil]. 
North is left and east at the bottom. Red hexagons show the zones of our IFS observations: (1) X-ray, (2) edge, (3) mini-bubble, 
and (4) bullet. 



ject. We have called all these regions : X-ray zone (1), edge (2), 
mini-bubble(3), and bullet (4) zones, as marked in Fig[T] 

We used three different gratings to obtain information for 
the most important emission lines in the optical spectral range. 
The R1200 grating used was centred in two wavelengths cover- 
ing the spectral range from 6220 to 6870 A and from 4320 to 
5060 A with a dispersion of 0.3A/pix, giving a spectral resolu- 
tion of -2.4 AFWHM at 6300A. A V1200 grating with a dis- 
persion of 0.35A/pix was used that covered the spectral range 
from 3660 to 4430 A with a spectral resolution of FWHM-2.7A 
at 4047A. This setup overlaps 100 A with the wavelength range 
of the R1200 blue grating. Finally, we used a low-resolution 
configuration with the V300 grating, covering the spectral range 
from 3660 to 7040 A with a dispersion of 1.67A/pix, giving a 
spectral resolution of -8.7 AFWHM at 5577 A. Table [2] shows 
the observational log for NGC 6888 in the four regions and 
three exposures were taken in each pointing. Moreover, in the 
X-ray zone with R1200 gratings, we adopted a dithering scheme 
with three positions having offsets measured from the centre 
of Aa = -1.35", AS = -0.78" for the second pointing and 
Aa = -1.35", A6 = +0.78" for the third. This pattern allowed 
us to fill up the gap between fibres. 



The weather was photometric throughout the observations 
with the typical seeing subarcsecond. Two spectrophotometric 
standard stars, Hz 44 and BD+284211, were observed for all 
the setups. The necessary bias frames, continuum, and arcs were 
acquired. 

2.2. Data reduction 

The data were reduced using the R3D software (Sanchez 2006) 
in combination with the Euro3D packages (Sanchez 2004) and 
IRAF Q The reduction consisted of the standard steps for 
fibre-fed IFS observations. 

A master bias frame for each night of observations was cre- 
ated by averaging the corresponding bias frames observed dur- 
ing the night and subtracted from all the images. The different 
exposures taken at the same position on the sky were combined 
to reject the cosmic rays using the package IMCOMBINE of 
IRAF. 



1 The Image Reduction and Analysis Facility IRAF is distributed by 
the National Optical Astronomy Observatories, which are operated by 
Association of Universities for Research in Astronomy, Inc., under co- 
operative agreement with the National Science Foundation. 
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Table 2. NGC 6888 PPAK observational log. 



Zone 


^UUrUlIlaLcS \ JZ.\J\J\J) 

(a,S) 


Grating 


Spectral range 
( A ) 


Exp. time 

(s) 


Airmass 


UdlC 


Comments 


X-ray zone (1) 


(20:12:39.8 , +38:25:34.5) 


V300 


3640-7040 


1350 


1.00 


July, 5 








R1200 


6200-6870 


1440 (4320 f ) 


1.02 


July, 1 


dithering (x3) 






R1200 


4320-5060 


1620 (4860 f ) 


1.03 


July, 2 


dithering (x3) 






VI 200 


3660-4430 


1620 (4860 f ) 


1.02 


July, 4 


dithering (x3) a 


Edge (2) 


(20:11:30.7 , +38:22:38.2) 


V300 


3640-7040 


1350 


1.01 


July, 5 








R1200 


6200-6870 


1620 


1.00 


July, 1 








R1200 


4320-5060 


1620 


1.00 


July, 2 




Mini-bubble (3) 


(20:12:29.0 ,+38:28:39.0) 


V300 


3640-7040 


1350 


1.00 


July, 6 




Bullet (4) 


(20:12:09.6 ,+38:20:14.9) 


V300 


3640-7040 


1350 


1.00 


July, 6 





' Total time including the three dithered exposures. 

" One of these dithering positions suffers from dome-lamp contamination. Only two of the three exposures are useful. 



Using continuum-illuminated exposures taken before the 
science exposures, we determined the location of spectra on 
the detector for each pixel along the dispersion axis. Then, 
each spectrum was extracted from the science and standard star 
frames, co-adding the flux within an aperture on the location of 
the spectral peak in the raw data using the tracing information, 
and storing in a 2D image called row-stacked-spectrum (RSS) 
(Sanchez 2004). Furthermore, we checked that the contamina- 
tion from flux coming from adjacent fibres using this aperture 
was negligible (Bershady et al. 2004 Sanchez 2006) . For a 
given aperture and FWHM ~ (0.5 X aperture) we found a 
level of cross-talk that is always <10%. This seems to be an 
acceptable compromise between maximizing the recovered flux 
and minimizing the cross -talk. 

Wavelength calibration was performed in a two-step proce- 
dure using arc calibration-lamp exposures (a neon lamp for the 
R1200 grating in the range from 6220 to 6870 A and He for the 
rest of the pointings). In the first step, distortion correction was 
carried out to avoid the curvature. Then, the wavelength coordi- 
nate system was determined by identifying the lines. Finally, the 
distortion and dispersion solutions were applied over the science 
data. The accuracy achieved was better than ~0.1A(rms) for the 
arc exposures in every spectral range. We performed a further 
wavelength "re-centring"using sky emission lines present in the 
observed wavelength ranges. 

Corrections to minimize the differences between 
fibre-to-fibre transmission throughput were also applied, 
creating a fibre-flat with an exposure of a continuum. All the 
science frames were divided by the obtained fibre-flat. 

Observations of spectrophotometric standard stars were used 
to perform the flux calibration. Only Hz 44 was observed all the 
nights, so we used it for our primary calibration. With the E3D 
software a spectrum for the calibration star was created each 
night by co-adding the central fibres of the standard star frames. 
Then, with IRAF tasks, a sensitivity function was determined. 
Subsequently, we checked that the flux-calibration for the two 
standard stars was self-consistent for all the nights observed and 
with the library values. Finally, flux calibration was applied to all 
the science frames. 

To achieve an accurate sky subtraction is an important issue 
in data reduction. Although PPAK has additional fibres for sam- 
pling the sky, our object is very extended, and all the sky fibres 
are located within an area containing signals from the nebula. 
Extra frames to sample the background near NGC 6888 were not 
available. For each night, we adopted different sky-subtraction 
schemes. For the night when we used the VI 200 grating, an ex- 
ternal sky-frame was observed. Since the sky and object expo- 



sures were taken within a few minutes of each other in the same 
night and under similar weather conditions, we decided to esti- 
mate a sky spectrum using this extra frame. To do so, the spectra 
of all individual fibres of the sky were combined with a mean in a 
single spectrum and a 2D spectrum was created with the value in 
each fibre of the combined spectrum. This sky was subsequently 
subtracted from every science spectrum. This method was used 
in observations with the high-resolution gratings (R1200 and 
VI 200), and no satisfactory sky frame was available for the night 
of V300 observations. 

Considering the short wavelength range and the airmass of 
the observations, and using Filippenko (1982), we estimated 



that the offsets due to the differential atmospheric refraction 
(DAR) were always smaller than one third of the diameter fibre. 
Correction for DAR was not necessary in our data. 



3. Two-dimensional study 

Once data had been reduced, we first carried out a basic 2D study 
to describe the observed morphology of the four pointings in 
several emission lines. Later, in zones with a dithering scheme, 
we performed a more detailed analysis. 



3.1. Morphology in the emission lines 

Using E3D software, we generated interpolated images of 
the four observed regions in three wavelength ranges in- 
cluding the lines: [OII]/U3726,3728A, [OIII]^5007A, and 
[NII]/U6548,6584A+H«. Figure||]displays these images. It can 
be seen two different patterns in all the regions, [Oil] and 
Hc+[NII] present similar behaviour, while [OIII] shows a dif- 
ferent structure. This can be understood as the signature of the 
ionization structure of the nebula. 

In the X-ray emitting region, most of the emission is lo- 
cated in the northwest and seems spread towards the northeast 
and centre. The high ionization emission (i.e. [OIII]) shows two 
main peaks, one in the northwest and the other in the southwest. 

The second zone (Edge) was chosen to be observed because 
is placed in a big bubble of [OIII] that appears surrounding all 
the nebula. It can be observed in Fig.|2j that in the [OIII]/15007A 
emission-line, a discontinuity separates the nebula into two dif- 
ferent regions, while it appears uniform in the low ionization 
emission-line images. 

The little bubble in the northern region of the nebula 
(Mini-bubble) is cleary seen in the [OIII] images, but appears 
very diffuse in Ha in Fig. [T] In the PPAK observations of this 
region different features can be distinguished, a bright zone to 
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(1) 



(2) 



(3) 



(4) 3 




Fig. 2. Images of the four pointings of NGC 6888. Columns are emission lines: from left to right [OII]/L?3726,3728A, [Offl],15007A, 
and [NII]/U6548,6584A+Hq'. The four different rows show the four regions observed; from top to bottom: (1) X-ray zone, (2) edge, 
(3) mini-bubble, and (4) bullet. In all the maps, north is down and east to the left. 



the south in all the wavelengths and a semi-circle structure only 
seen in the [OIII] image. 



Finally, the last row of Fig. ^corresponds to a darker region 
placed near the centre of NGC 6888 that could be a sort of 
"bullet" crossing the nebula. This region appears quite similar 
in the three wavelength maps. The dark zone in the north is very 
faint in all the emission lines, though again low ionization lines 
show slight differences from [OIII]. 



3.2. Data-cube analysis of the X-ray zone 

The bi-dimensional structure of the three dithering observations 
was studied in detail by means of the creation of data cubes. 
We performed this kind of analysis only for the X-ray zone as 
observed with R1200 gratings in two wavelength ranges because 
for the other pointings dithering was not available. In the next 
sections (3.2.1 3.2.2 and 3.2.3 i, all the results refer to the X-ray 
zone. 

As explained above, the 331 science fibres lead the light 
from different positions of the target to the pseudo-slit. When 
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dithering exposures (to fill holes) were performed, it is impor- 
tant to combine them when knowing the overlapping areas and 
registering the spectra to their positions in the sky. With this 
aim, we used the R3D software to construct a grid with a pixel 
scale about one third of the original fibres (~1 arcsec), obtaining 
easier-to-use data cubes with two spatial and one spectral 
dimension. Therefore, two cubes of the X-ray region in the 
wavelength ranges 6220-6870 A and 4320-5060 A were built. 

One of our objectives was to obtain maps from the datacubes 
of the main emission lines, which could then be analysed to 
describe the 2D ionization structure and physical properties of 
the X-ray emitting region of NGC 6888. This procedure was 
performed using our own routine, which fits each line with a 
Gaussian function. We introduce four parameters as a reference 
to the fitting (the central wavelength, the FWHM, the line wave- 
length range, and the corresponding continuum range) and maps 
are returned with useful information in each pixel for each se- 
lected line (i.e. flux, FWHM, central wavelength, error of flux, 
continuum, etc.). 

As a first step, we fit all the lines by a single Gaussian for 
the two data cubes obtained. For the grating placed in the red 
range (from 6200 to 6870 A), seven lines were fit: [Slll]A 63 12A, 
Ha, [NII]X16548,6584A, HeL16678A, and [SII]^6717,6731A. 
In the other spectral range (from 4320 to 5060 A), another 
seven lines were fit: Hy, UB, [OIII]/L14363,4959,5007A, and 
HeI/l/l4471,5015A. In addition, to prevent contamination by low 
signal-to-noise (S/N) data, we masked all pixels with an S/N 
lower than 5 out. Under these conditions, maps were used in the 
analysis that we present in the next sections. 

Two main problems appeared in the maps creation. First, 
with those fainter lines placed in zones with high noise, the fits of 
the automatic routine were not good in all the spaxels, and some 
adjacent signal was added to the fit, giving a wrong measure- 
ment of the line flux. It happened with three lines: [SIII]/i6312A, 
HeL14471A, and HeLt5015A. With the obtained maps for these 
lines, masks were created in order to not take bad pixels into 
account in the analysis. The second problem that we found was 
fitting the [OIII]/!4363A. This line is located very close to the 
strong Hgl/14358A night sky line, in particular, both lines are 
overlapped in our object, and we cannot resolved them, so we 
could not obtain workable maps of [OIII]/l4363A for this neb- 
ula. Without the measure of this line, it was not possible to ob- 
tain the electron temperature map. Furthermore, no other auroral 
line, like [NII]/15755A, was measured in the cubes with a high 
S/N (but see Sect. lO). 



3.2.1. Extinction and electron density maps 

The reddening coefficient c(HB) was derived from the Ha/HB 
line ratio. We used an intrinsic Balmer emission line ratio of 
Ha I Hp = 2.81 obtained from the public software of Storey& 
Hummer (1995) assuming Case B recombination with an elec- 
tron temperature of T e ~ 10 4 K and an electron density of 
« e ~ 10 2 CUT 3 . 

Although other Balmer lines were detected in our wave- 
length spectral range, Ha and HB were the lines with the highest 
S/N. We performed the estimation of c(Hy8) including Hy. The 
average value obtained was similar but the errors increased, and 
the S/N in the Hy flux was lower than 5 in 10% of the spax- 
els, so finally we only used Ha and HB. However, we checked 
that in the spaxels with good S/N in Hy both derivations were 





Fig. 3. Reddening coefficient map, c(HB), derived from the 
Ha/HB line ratio. Orientation as in Fig. [2] 



consistent. The derived c(HB) map is shown in Fig. [3] where a 
non-uniform, rather filamentary structure can be clearly distin- 
guished, with values ranging from c(H/3) = 0.3 to 1.00, with a 
mean val ue of 0.57. Adopting E (B - V) = 0.692 x c(H/3) and 
using the Car delli et al.| ( T989] > extinction law with R v = 3.1, 
E(B-V) was computed with values in the range 0.21-0.88. In 
general terms, the structure presented in this map is in good 
agreement with values reported in the literature. Using the flux 
density at 21cm and the flux in Ha, Wendker et al. ( 1975) esti- 
mated an absorption over the whole nebula of Ay ~ 2.05 mag 
(E(B - V) = A v /3.1 = 0.66). |Esteban & Vilchez| ( fl992| ) ob- 
tained values of the c(Hy6) reddening coefficient of 0.28 for the 
blueshifted component and 0.50 for the redshifted one, corre- 
spond to E(B-V) of 0.19 and 0.34, respectively. Some of these 
values are slightly lower than ours, but their slit positions do no 
spatially coincide with our observed zones. 

Finally, all the emission-line maps were 
reddening-corrected using the derived c{HB) map. Observing 
their appearance we checked the morphology explained in Sect. 
3~T]but now with more lines: Ha, HB, [SII], [Nil], and He I with 
similar patterns as examples of low ionization and the doublet 
of [OIII] as high ionization. 



The electron density (n e ) map of the X-ray region was 
produced from the [SII]/l/l6717/6731 ratio using the package 
TEMDEN of IRAF based in a five-level statistical equilibrium 
model ( |De Robertis et al.|1987|[Shaw & Dufour|1995) ). The spa- 
tial distribution of the electron density can be seen in Fig. [4] The 
map shows a rather filamentary aspect that is more or less uni- 
form with a mean value of 178 crrT 3 . The main significant fea- 
ture in this spatial distribution is the peak density (~400 cm -3 ) 
coincident with the maximum flux emission of all the species 
detected. Our highest n e values agree with the density obtained 
by |Kwitter| ( [T98T) and jEsteban & Vilchezl ( [T992 i. 



3.2.2. Diagnostic diagrams 

The excitation mechanisms and the ionization structure have 
been studied using diagnostic diagrams based on the ratio of 



strong emission lines: the classic BPT analysis (Bal dwin et al.| 
1981 ). Using the data from the 2D study of the X-ray zone, we 
generated maps for the necessary line ratios (with S/N > 5) and 
plotted all the spaxels in the three diagrams displayed in Fig. [5] 
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Fig. 4. Electron density n e map derived from the 
[SII]/i/i67 17/6731 line ratio. The n e distribution ranges from 
~50 cirT 3 to a maximum peak of 410 cm -3 , with a mean of 178 
cm 4 . Orientation as in Fig. [2] 



The first diagram (Figga) shows [NITJ/U6584/Ha vs. 
[SII]/1/16731/Hq' (NvsS) and presents a bivaluate behaviour. We 
performed a reference fit to the points, paying attention to which 
spaxels are inside/outside the nominal limits of +3<x from this 
fit. It can be seen that the relation between these two line ra- 
tios shows a significant scatter, especially for the points below 
the -3cr line. We notice that spaxels above the +3cr line show 
a linear-like correlation as could be expected from the ioniza- 
tion structure given the ionization degree of these ions. To un- 
derstand these behaviours, we located the points of the diagram 
in the map differentiating between points inside, above and be- 
low the 3cr limits of the global fit, to identify their spatial loca- 
tions (see Fig j6]». We found that points with a similar tendency in 
the BPT appear grouped in the PPAK map. The second diagram 
(Fig.[5]b) presents [OIII]/U5OO7/H0 vs [SHJ/U673 1/Hor (OvsS), 
it shows a larger scatter with two tendencies. To understand how 
the zones defined above in Fig j6] behave in this diagram, we 
identify pixels belonging to these zones with the same colours 
as in Fig. |5]a. Points with large scatter in NvsS (blue crosses) 
show scatter also in this diagnostic diagram, and points with a 
correlation in NvsS (red crosses) appear here an ti -correlated. 
In the last diagram (Fig. |5]c) we can see the relation between 
[Oin],U5007/H/? and [NII]X*6584/Har (OvsN), where we find a 
strong anti-correlation. Once again, we study the behaviour of 
the three regions defined in Fig|6]in this diagram with the same 
colour code. In this case, the points appear anti-correlated for 
all the zones. 

Taking the analysis and results explained above into account, 
we have defined three spatial zones in this pointing: Zone A to 
the northwestern part of this pointing, corresponding to all the 
pixels above +3cr in NvsS. We notice that in all the diagrams 
these spaxels show a normal and expected tendency. In addition, 
in this area the maximum emission in low excitation species is 
found. Zone B is located to the southwest and corresponds to 
points below -3cr in NvsS. In this case, the relations between 
the line ratios in the two first diagrams were rather anomalous 
and had a large dispersion. In the map this zone corresponds 
to the peak of the [ODTJ/15007A emission. Finally, we defined 
Zone C by the points inside ±3<x limits of the fit in the diagram 
NvsS, which are located over the rest of the field and present a 
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Fig. 5. Diagnostic diagrams with data from the 2D study of the 
X-ray pointing. From the top to the bottom: (a) [NII]/l/16584/Ha 
vs [SII]/Li6731/Ha, (b) [OIII]i/l5007/HyS vs [SII]/Ll6731/Ha<, 
and (c) [OIII]ii5007/H/3 vs [NII],U6584/Ha. All the spaxels 
of the intensity maps are represented with crosses. Black lines 
in the first diagram are the fit performed and the +3cr limits. 
Colours help us to locate points spatially in Fig. [6] red corre- 
sponds to Zone A, blue is Zone B, and green represents Zone C 
(see text for details). 



behaviour that is a mixture between the two other zones. 



3.2.3. Kinematics 

The radial velocity of the gas was derived using the central wave- 
length of the Gaussian fit in every pixel. All the measured radial 
velocities were translated into the local standard of rest (LSR) 
and corrected for the Earth's motions using IRAF. In addition, a 
zero-velocity map was created using a "sky" line, and this map 
gives us realistic information about the error in the wavelength 
calibration and the accuracy in the velocity that can be achieved. 
We concluded that differences in radial velocity below 8 km s _1 
cannot be covered with our data. Figure [7] shows the velocity 
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Fig. 6. PPAK FoV of the X-ray pointing with the different ion- 
ization zones defined by taking the tendencies in the BPT dia- 
gram of Fig. [5]a into account. Red corresponds to Zone A, blue 
is Zone B, and green represents Zone C. Orientation as in Fig. [2] 




Fig. 7. Radial velocity field derived for Ha. Velocity ranges from 
60 to 1 10 km s _1 with a mean value ~84 km s _1 . Orientation as 
in Fig. [2] 



field derived from the Ha emission line. The map has a mean 
value ~84 km s _1 , but it is not spatially homogeneous, because 
the southwest zone moves faster, reaching 110 km s _1 , while in 
the northwest zone the velocity is lower than 60 km s _1 

Taking advantage of the large number of pixels in our maps, 
we have produced statistical frequency distributions of the radial 
velocity. When we represent velocities from all the pixels, the 
histograms present a Gaussian distribution as expected. But, if 
pixels from Zones A and B are represented separately, a bimodal 
distribution is found in some relevant emission-lines like 
[NII]/16584A. Figure [8] shows two histograms generated with a 
binning of 3 km s _I for two emission lines differentiating the 
two defined zones. The [NII]/16584A line presents two patterns: 
while Zone A shows a distribution with a peak redshifted by 
10 km s , Zone B shows a bimodal distribution indicating a 
shell-like behaviour with pixels showing gas that is moving 
away and other pixels coming towards us. On the other hand, in 
the second histogram, the [OIII]/i5007A line presents a single 
peak distribution in both zones and the velocity difference 
between the two peaks is ~15 km s _1 . 

With this kinematic statistical study we support the idea that 
NGC 6888 is a non-uniform nebula with shells and filaments 




Relative radial velocity (km/s) 



Fig. 8. Statistical frequency distributions of the radial velocity of 
[NIIJ/16584A and [OIII],t5007A. Solid red lines indicate pixels 
from Zone A and dashed blue lines from Zone B. Velocities are 
given in units of km s taking as zero the mean velocity mea- 
sured for the emission line. 



moving with different velocities. Errors in the wavelength cali- 
bration and the instrument resolution of the observations prevent 
us from carrying out a more exhaustive analysis of the kinemat- 
ics of the nebula in this work. 



4. One-dimensional study 

A different analysis of the IFS data can be carried out by 
combining the observed spectra in several fibres to produce 
integrated spectra in different regions of NGC 6888. With 
this method we can use data without dithering exposures to 
obtain information over a broader spectral range and derive 
some physical properties that could not be obtained using a 
2D approach. To decide which fibres are to be combined, we 
examined the appearance of each zone in emission lines at 
several wavelengths using E3D as shown in Fig.|9](see also Fig. 

- X-ray zone: two integrated spectra were created in this point- 
ing. The first spectrum (XI) corresponds to spaxels where the 
[NII]/15755A line is visible. We have chosen this criterion be- 
cause we were looking for emission lines useful for deriving the 
electron temperature. In addition, spaxels where [NII]/15755A is 
bright enough coincide with the brightest low-excitation lines. 
To define the second region (X2), we selected spaxels placed 
in the area where we had a peak in the high ionization line 
emission ([OIII]). Moreover, XI and X2 are in Zones A and B, 
respectively. 
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Fig. 9. Displays using E3D of the four areas observed with 
PPAK. Boxes show regions where the integrated spectra were 
created: (1) X-ray zone, (2) Edge, (3) Mini-bubble, and (4) 
Bullet. All the images are RGB -combinations of three wave- 
lengths: red corresponds to the Ha+[NII], green to [OIII], and 
blue to [Oil]. Arbitrary scales. Orientation as in Fig. [2] 



- Edge zone: we were interested in checking the different condi- 
tions and parameters on both sides of the nebular discontinuity 
observed in [ODTJ/15007A, so the first integrated spectrum was 
performed in the brighter zone to the left (El) and the second 
one in the fainter part (E2). 

- Mini-bubble: we created integrated spectra over three of the 
features seen. Firstly, we chose fibres emitting in [NII]/15755A 
(MB1). Secondly, a bright zone in [OHI]/i5007A but without 
[NII]/15755A emission (MB2). Finally, we defined the zone 
with low level of emission at all the wavelengths (MB3). 

- Bullet zone: to check the nature of the bullet candidate, we 
created one integrated spectrum over the dark zone (Bl). In 
addition, as a representative example of the emission of the 
inner area of the nebula, we extracted a spectrum (B2) in this 
brighter region but outside of the bullet. 

Using these definitions, several sets of fibres were selected 
and combined for each grating in all the zones to obtain their 
integrated spectra, which were analysed to derive their spectro- 
scopic properties in the whole spectral range of the observations. 
For observations where a dithering scheme was adopted we com- 
bined the fibres of all the exposures. At the end, we obtained 19 
ID spectra considering the different gratings. Figure 10 shows 
an example of integrated spectra in zone XI. 



4. 1 . Integrated properties 

The emission-line fluxes were measured using the SPLOT task 
in IRAF. Lines show a degree of asymmetry in their profiles, so 
we did not perform Gaussian fits The fluxes were derived mea- 
suring areas below the lines, however, partially blended lines, 
such as Ha and [NII]/U6584,6584A in the V300 grating, were 
deblended with two or three Gaussian profiles to measure the 
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Fig. 10. Integrated spectra corresponding to the low ionization 
area of the X-ray zone (XI). Top: the whole spectrum in the 
V300 grating. Middle: a zoom centred on the [NITJ/15755A and 
Hel/15876A emission lines. Bottom: the R1200 red set up. 

line fluxes separately. The statistical errors associated with the 
observed fluxes were calculated using the expression 



o-j = cr c N l,2 [l +EW/(NA)] 



1/2 



(1) 



where crj represents the error in the observed line flux, N is 
the number of pixels used to measure the line, EW the line 
equivalent width, cr c the standard deviation of the continuum 
in a box near the line, and A represents the dispersion in A/pix 
flPerez-Montero & Diaz|2003) . 

Table [3] lists the emission lines measured in every zone la- 
belled with its standard identification. The estimated fluxes and 
errors were normalized to F(HB) = 100. Line intensities were 
reddening-corrected using the same procedure as described in 
Sect. |3.2.1| In this case, other HI Balmer lines with S /N > 5 
were measured in the integrated spectra and therefore used to 
derive the reddening coefficient. In all the spectra, Ha, HB, Hy, 
and H(5 were measured with enough S/N and used to estimate 
c(HB), with the exception of the pointing Bl in which Hy and 
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H6 were very faint. The reddening coefficients derived in zones 
XI and X2 agree with the values obtained in the 2D-study in 
the X-ray pointing using only Ha and H/3. The values obtained 
for c(H/J) are presented in the last row of Table [3] Errors in the 
emission-line intensities were derived by propagating the ob- 
servational errors in the fluxes and the reddening constant un- 
certainties. 

In regions XI, X2, El, and E2, some lines were measured 
more than once since they appear in exposures with different 
gratings. In these cases, we calculated the corresponding 
intensity as the mean weighted by the error in each zone. 

To obtain the physical conditions of the gas, we performed 
an iterative process until achieving agreement between electron 
temperature and electron density. Table |4] shows the values for 
T e and n e derived for each zone. 

With the [SITJ/l/167 17/6731 line ratio, electron densities (n e ) 
were obtained using the package TEMDEN of IRAF. We derived 
differences in the electron density of up to 350 cirT 3 between 
the studied zones. The results from the X-ray zones are coherent 
with the values obtained in the bi-dimensional analysis where 
densities were higher in the northwest of our field of view that 
corresponds to the integrated spectrum XI. On the other hand, 
values obtained on the edge of the nebula (El and E2) and in 
the less bright zones of the mini-bubble (MB2 and MB3) are 
lower than the densities derived in other regions (always n e < 
100 cm" 3 ). 

Electron temperature, T e , can be derived using the line ratios 



the forbidden-to-hydrogen emission-line ratios using the func- 
tional forms given by Hagel e et al.| (2008), which are based on 
the package IONIC of IRAF. To determine the singly ionized 
helium abundance, we used the equations proposed by Olive 



7?03 



and 



R N 2 - 



/([OIIIM4959) + 7([OIII]/15007) 
7([OIH],14363) 



7([NII]/i6548) + 7([NII]^6584) 
7([NII]^5755) 



(2) 



(3) 



The [OIII]i4363A line was blended with the "sky" line 
HgL14358A, and therefore we cannot compute T e ([OIII]) di- 
rectly. In contrast, the [NII]/15755A auroral line was measured 
with enough confidence in zones XI, MB1, and B2, and di- 
rect estimates of T e ([NII]) were obtained from the Rn2 ratio. In 
the rest of zones (X2, El, E2, MB2, MB3, and Bl), for which 
this line was not bright enough to be measured with confidence, 
we resorted to the empirical flux-flux relations from Pilyugin 
(2005) that relate the faint auroral lines with other strong neb- 
ular lines. To find out what was the best relation to our partic- 
ular case, we used several fits proposed by Pilyugin over the 
zone XI, for which Rn2 is well measured, and selected which 
relation represented our results better. We chose the relation 
Gnu = f(Ri,P) from |Pilyugin| ( |2007| l to estimate the elec- 
tron temperature T e ([NII]) and checked that this relation gives 
consistent results with our measurements of the line in the XI 
zone. Later, from T e ([NII]), we derived T e ([OIII]) and T e ([OII]) 
in each zone by using the relations based on models proposed 
by Perez-Montero & Contini ( 2009 \ and Perez-Montero & Diaz 
(2003|), respectively. Derived temperatures range from ~7700 K 
to ~ 10 200 K over the nebula, which is consistent with the pre- 
vious studies of |Kwitter| ( | 1 98 1| > and |Esteban & Vilchez| ( |1992| >. 

4.2. Chemical abundances 

Considering the electron density and temperatures estimated 
for each ionization zone, we obtained ionic abundances from 



& Skillman] ( |2004| l. We defined two temperatures: T e ([NII]) as 
temperature representative for the low ionization ions, S + , N + , 
and + , and T f ([OIII]) for the high ionization species, 2+ and 
Ne 2+ , and for deriving the abundance of He + . All the ionic abun- 
dances are presented in Table[4]with their corresponding errors. 

Furthermore, in Table |4] we show the total abundances de- 
rived as follows. The O/H has been obtained by directly adding 
the two ionic abundances (0 + and 2+ ). The same procedure 
was used for He/H but only with He + because He 2+ is not seen in 
our observations, the value chosen of Hel was a mean weighted 
by the errors of the three Hel observed lines (HeL14471A, 
HeLi5876A, and HeL16678A). The N/O can be approximated to 
N + /0 + . Finally, to derive the total Ne/H abundance, ionization 
correction factors (ICF) are required to take the unseen ioniza- 
tion stages into account, we used the ICF expression for Ne from 
Perez-Montero et al. ( 2007[ l. Taking Izotov et al. ( 2007[ ) as refer- 
ence, we estimated ICF(He) ~ 1 for all the zones, so correction 
for the unseen neutral helium was not necessary. 

Mainly, three results in the total abundances are interesting. 
First, we see that nearly all the zones present an oxygen abun- 
dance slightly below the solar value, taking the solar abundance 
from |Asplund et aL] ( |2009| ) as reference . This effect was reported 
above for the areas observed by |Kwitter ( 198 1| > and |Esteban &] 
|Vilchez| (fr992 ). On the other hand, these authors derived an over- 
abundance of nitrogen. For some of our zones (El, E2, MB2, 
MB3, and Bl), we derived a nitrogen abundance similar to the 
solar value, however in other zones the N/H appears enhanced 
by a factor of 6, or even 8 in the XI zone. Helium presents an 
enrichment in most of the integrated zones too. (MB3, E2, and 
Bl do not show the Hel line.) One of the most striking results is 
the important differences in the N/O between the different zones, 
mainly due to the large N/H variation and not to differences in 
O/H. This could be understood by assuming that we are seeing 
structures of the nebula formed in different stages of the evo- 
lution of the WR star, each one with its own ionization degree, 
physical properties, and chemical patterns. This idea is analysed 
in Sect. [5] Furthermore, we calculated the expected values in the 
ISM at the galactocentric radius of NGC 6888, considering the 



galactic chemical abundance radial gradients. Taking Balser et 
[aLl( [20rT] ) as reference, we derived 12 + log (O/H ) = 8.6 + 0.05 



and log (N/O) = -1.26 + 0.26 from |Carigi eUriTlpOOB) . Our in- 
ferred abundances for E2, MB3, and B 1 are consistent with these 
expected values (within the errors). 

4.3. Photoionization models 



To investigate and consolidate both ionization parameter (U) 
and chemical abundances in the different pointings of the neb- 
ula, we performed a grid of photoionization models made with 
the code Cloudy v.8.0 (Ferl and et al.||1998[). These mod els are 
based on the grid described by Perez-Montero & Contini (2009 ), 
using one single stellar atmosphere of WR, in particular, a 
WN star model ( |Smith e t al. 2002). Our grids assume an ef- 
fective temperature of the star of 50 000 K ( |Rosa & Mathis 
19901 |Esteban et aI1[T993) , with a meta llicity of Z = 0.008 
Taking the solar abundances reported by Aspl und et al.| ([2009 ) 
(12 + log (O/H) = 8.69, Z = 0.013) as a reference, the assumed 
metallicity implies an oxygen abundance 12+log (O/H)) = 8.47, 
which is the closest value to the O/H abundances derived in 
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Fig. 11. Diagnostic diagrams containing the results 
of our photoionization models and observed data, (a) 
[OIII],U5007/H/? vs [Nn]/U6584/Hor and (b) [OIII],M5007/H/3 
vs [NII],16584/[OII]X13726,3729. In all the diagrams, the dotted 
lines connect points with the same ionization parameter. The 
dashed lines connect curves of equal log(N/0). The respective 
values are indicated at the beginning of curves. For all the 
models we assumed 12 + log(0/H)) = 8.47. Points are values 
of integrated spectra placed according to their intensities ratios 
(named as in Tableland Pig |9j. 



this work for NGC 6888. These models vary the ionization pa- 
rameter and the nitrogen-to-oxygen ratio as free parameters, 
which are very relevant for our purpose. We considered values 
log(f/) = -2.5, -3.0, and -3.5, and log(N/0) from -1.5 to 0.0 
(in steps of 0.25). 

We compared selected emission-line ratios of the in- 
tegrated spectra with models in two diagnostic diagrams: 
[OIII],U5007/H/? vs [Nn]/U6584/Ha and [OIII]/U5007/H/3 vs 



[NII],16584/[OII]Xl3726,3729 (Dopita et aDgOOO]). As can be 



seen in Fig. 1 1 for each zone, the values of log (N/O) obtained 
with the comparison always agree with abundances calculated 
with our measures. For each integrated spectra, the ionization pa- 
rameter is similar in the two diagrams: zones X2, El, and MB2 
have the greater U; zones XI and B2 have intermediate values; 
and the paremeters are the lowest in zones E2, MB1, MB3, and 
B1.0 



2 Models with an effective temperature of 45 000 K and WM-Basic 
stellar atmosphere l Pauldrach et al. 200 lj l also gave satisfactory results. 




Fig. 12. Sketch showing the structure of NGC 6888 with two 
main shells: one inner (IS) and the other outer (OS). 



5. Discussion 

To develop a scheme of the internal structure of NGC 
we put together all the information derived from the 2D maps 
and the nine integrated spectra. Furthermore, we used the 
narrow-band images from the INT-WFC shown in Sect. [T]to 
have a global overview of the nebula and to locate regions of 
different observed conditions. The proposed scheme consists of 
a shell structure (close to "onion-like") but allowing irregular 
and/or broken shells. We found, at least, two shells with dif- 
ferent physical, kinematical, and chemical properties: an outer 
shell (OS) and an inner shell (IS). In addition, a more external 
and fainter structure is seen surrounding the whole nebula like a 
thin skin; we call it surrounding skin (SS). A basic sketch of this 
scheme is displayed in Fig [12] 

Our interpretation of the structure of NGC 6888 is supported 
by comparison to models of WR ring evolution. In particular, 
there is very good agreement with the simulation by |Garcia- 
Segu ra et al.| JT996} (hereinafter GS96) for the interaction of 
a 35 Mq star with the circumstellar medium over the entire 
stellar lifetime in the case of the slow RSG wind (~15 km s _1 ). 
Their simulations start at the MS stage, when the stellar winds 
sweep up the ISM, forming a hot bubble surrounded by a 
thin shell. Later, in the red supergiant (RSG) phase, the star is 
characterized by a slow and very dense wind that shocks with 
the MS bubble, creating a shell (RSG shell). Finally, a faster 
wind from the WR phase sweeps up the RSG wind material 
into a new shell; however, the propagation of this WR shell is 
faster, and it soon arrives at the RSG shell. After the collision, 
clumps and dense material of the WR shell are propelled by 
their high inertia, crossing the RSG shell as bullets forming 
blow-outs. The final stage in the GS96 simulations (Fig. 7 in 
that paper) shows that the bubble formed in the main sequence 
is, eventually, broken out as consequence of the RSG-shell 
shocked by the WR shell and its expansion. 



5.1. The inner shell (IS) 

In the proposed scheme, the central part of NGC 6888 is an 
ellipsoidal shell surrounding the WR star with a high surface 

T). This shape 



brightness, as can be seen in the INT image (Fig 
agrees with the scheme proposed previously by |Lozinskaya et 
[ir]([T988). Based on interferometric observations, they derived a 
spatial model in which NGC 6888 is a prolate ellipsoid of revo- 
lution (a/b ~ 0.5) whose major axis is inclined 20° - 40° to the 
projection plane. 
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In our data, the information about this shell was obtained 
from the integrated spectra XI, X2, MB1 and B2. From the 
analysis of the properties of these four zones, we propose that 
this shell is not a complete shell, but a broken structure, with 
holes, filaments and bubbles. One evidence pointing to this are 
the variations of the electron density always >100 cirT 3 ranging 
from 106 to 360 cirT 3 along the defined IS. Also, differences in 
the electron temperature and in the chemical composition show 
inhomogeneities between different zones of this shell. 

One of the more important results for the IS is the strong 
overabundance in N/H by a factor between 5 and 9 with respect 
to the solar abundance and by a factor up to ~11 with respect 
to other regions of the nebula. The enriched N/H found varies 
from 12 + log (N/H) = 8.29 to 8.72 in the integrated spectra in- 
cluded in this shell, supporting the idea of a patchy and incom- 
plete structure. On the other hand, the O/H abundance is found 
to be slightly deficient, suggesting that some oxygen could have 
been processed into N mainly via the ON cycle. In the case of he- 
lium, we obtained a value of up to He/H=0.18 in this shell, indi- 
cating a strong overabundance of this element when we compare 
it with the solar value He/H=0.09 ( |Asplund et al.|2009| >. These 
facts give us clues to the stage in the star evolution in which this 
shell originated. 

We have additional information about the IS thanks to 
the observations with IFS because the X-ray pointing (with 
dithering and full datacubes) is over this region. In Sect |3.2.2 
we defined two zones, Zones A and B, associated with the 
emission of the low and high ionization ions, respectively, and 
showing very different trends in the diagnostic diagrams. Later, 
we found two behaviours in the radial velocities associated to 
these zones. In particular, the statistical frequency distributions 
of velocity (Fig. [8} are consistent with a broken shell scheme: 
the [NII]/16584A velocity histogram in Zone A suggests only 
one redshifted shell going away from us, while in Zone B the 
bimodal velocity distribution found can be understood if the 
shell is moving away in some places and in others is coming 
towards us. Rarely are the two behaviours observed simultane- 
ously, because we do not see a central and dominant velocity 
component in the histograms. This kinematical study supports 
the idea of an intrinsically inhomogeneous shell in expansion, 
more obvious in the nitrogen than in the oxygen emission. 
Figure [13] shows a sketch of how the scheme of broken-shell 
works and what we would see for each zone. 



Although in our work we cannot consider the 2D study of the 
whole nebula, including all the spatial information, we propose 
that the most internal shell of NGC 6888 should present a struc- 
ture with inhomogeneities in the physical properties and abun- 
dances. Differences in the chemical composition could be under- 
stood by attending to the phase in the evolution of the central star 
when the shell was formed; in particular and in comparison with 
GS96 simulations, we argue that the IS material corresponds to 
the red supergiant (RSG) and WR shells shocked expanding in 
a low-density medium. This would explain its fragmented and 
filamentary structure and the high emission of Ha. Also, the un- 
usual abundances found indicate that this shell is mostly com- 
posed of stellar material ejected from the central star when it 
became a WR star and later probably mixed with material from 
the RSG stage. This idea is also supported by the recent simu- 
lation performed by Toala & Arthur (2011). They find that the 
excess of N/O in NGC 6888 can be explained and reproduced by 
a stellar evolution model with rotation (Meynet & Maeder 2003| > 
for an initial stellar mass of 37 Mq. 
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Fig. 13. Sketch showing the proposed broken inner shell. When 
we look at Zone A, we only receive light from the furthest shell 
that is moving away from us (the histogram shows a redshifted 
peak); on the other hand, when we look at Zone B, we see gas 
belonging to both shells (blueshifted or redshifted) but not si- 
multaneously or necessarily in the same spaxels. 



5.1 .1 . Can we detect the footprint of shocks? 

The previous observations of X ray emission detected in this re- 
gion and the theoretical simulations (GS96) lead us to consider 
that on the north edge of the major axis of this IS, a collision 
between the WR and RSG shells happened during the star's life- 
time. To check that there are spectral signatures of shocks, we 
resorted to the models presented by |Allen et al. (20081. Their 
code, 



MAPPINGS III, gives an extensive library of radiative 
shock models covering a wide range of abundances and veloci- 
ties. First, we had the maps of the emission lines ratios created 
in our 2D study; then, we generated a grid of shock models with 
different parameters using libraries from |Allen et al.| ( |2008[ ); fi- 
nally, we overplotted all the spaxels from the maps on the the- 
oretical grids in the diagnostic diagrams to find out the consis- 
tency with models. 

In a first attempt, we plotted all the pixels of our 
maps in two diagnostic diagrams ([OHI]/i/15007/Hj8 
vs [NII]/Ll6548,6584/Ha and [OJll]AA5Q07 /HB vs 
[SII]/L16717,6731/Ha<), but we failed to find a model that 
represented all the observational data in both diagrams well. 
Later, we studied only the region with a striking behaviour in 
Figs. [5 and[6](called Zone B). When we plotted only the spaxels 
that belong to this region, we found that they can be fit with one 
of the grids of models. Figure 14 represents the two diagnostic 
diagrams showing the pixels from the Zone B, together with the 
grid of shock models. 

We found that our data of Zone B can be represented well 
by a model in which shock and precursor are involved assuming 
n=1000 cirT 3 and with twice solar abundances. In both diagrams 
we find that the models that fit better have shock velocities from 
250 to 400 km s _1 . The agreement of our data with the predic- 
tions of shock models with twice solar abundances can be un- 
derstood by attending to the strong overabundance of nitrogen. 
Nevertheless, we should bear in mind that in our ID study the 
line ratios of this region (integrated spectra X2) are also consis- 
tent with the tailored photoionization model presented in Sect. 
14.31 An alternative scenario can be considered that is similar to 
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Fig. 14. Diagnostic diagrams: [Om]/U5OO7/H0 vs 
[NII]^^6548,6584/Ha top, and [Oin]X15007/H/? vs 
[SII]/l/i67 17,673 1/Ha bottom. Grids represent the predic- 
tion from the shock models using libraries from | Allen et al.| 
(2008). Numbers over the lines indicate the value of the shock 
velocity represented by red lines, while the blue lines are 
models with equal magnetic field. Crosses are the pixels from 
NGC 6888 maps belonging to the Zone B. Only data from 
Zone B are represented well in both diagrams by a model 
with shock and precursor, velocities from 250 to 400 km s _1 , 
n = 1000 cm 3 , and twice solar abundance. 



the ionization of "pillars of creations" structures (Ercolano et al. 
201 l| l. In this scenario, altered line ratios could be produced as 
a consequence of the partial observation of ionized to neutral 
surfaces in the nebula. 

As a conclusion, considering the shock models, the X-ray 
emission and the stellar evolution simulations, we cannot rule 



out the presence of a shock, so if this event has happened, the 
analysis of our data indicates that the most likely place to find it 
is Zone B of the X-ray pointing. It would be interesting to check 
this with new observations of the other side of the major axis 
(southwest of NGC 6888). 



5.2. The outer shell (OS) 

Studying the integrated spectra El and MB2, we found several 
differences in the physical properties in comparison with other 
zones that lead us to think that they belong to a different shell. 
In our nebular scheme, we propose an intermediate shell placed 
between the SS and the IS. The morphological study made in 
the bi-dimensional analysis (Sect. |3.1| l supports the existence 
of different shells or spatial inhomogeneities in the structure. In 
particular, the images in the second row of Fig. [2] correspond- 
ing to the Edge pointing, show a clear discontinuity in the spa- 
tial distribution of the emission of the [OETJ/15007A line that 
is not visible in Ha nor [OH]/i3727A. The integrated spectrum 
El was created from the emission to the right of this disconti- 
nuity. Something similar happens in the mini-bubble pointing 
(third row in Fig. [2]) where a well-shaped semicircular struc- 
ture (where the MB2 spectrum was extracted) is seen only in the 
[OIII] emission line. 

Taking the previous observations of NGC 6888 into account 
( |Gruendl et al.|2000| [Moore et al.|2000| l, it is very likely that the 
OS has an apparent spherical shape, at least in the [OIII]/15007A 
emission line. The INT-WFC image (Fig. [TJ agrees with this 
morphology: it can be seen that the oxygen distribution is not as 
elliptical as the other elements, showing some extended emission 
beyond the end of the minor axis of the IS. 

This bubble is surrounding the inner shell and looks domi- 
nant in double ionized oxygen. Intensities measured of the emis- 
sion lines are weaker here than in the inner shell; nevertheless, 
the [OIII]/15007A line has the maximum flux in comparison with 
the rest of the integrated spectra. Here we estimate the highest 
electron temperatures of the whole sampled nebula (>10 000 K). 
The expansion of the shell could explain the electron density 
found here, always lower than 100 cm 4 . The chemical anal- 
ysis shows that the He/H is enriched with respect to the solar 
values, while N/H and O/H appear not to be enhanced, or they 
present a slight deficiency. It is interesting to notice that the O/H 
abundance is similar to the one of the IS, but N/H is lower by 
a factor ~8. It was not possible to estimate the neon abundance 
here because these lines were not bright enough in the integrated 
spectra. 

We inferred the gas distribution in this shell studying data 
from the maps in the X-ray zone, because, although in our 
scheme the X-ray pointing does not target the OS, it gives 
us other interesting information. Remembering the statistical 
frequency velocity distributions in Fig. [8] it can be seen that 
the overall range in radial velocity covered by [OIII]/15007A 
and [NII]/16584A emissions is similar (~80 km s ), but the 
[OIII]/15007A velocity distribution shows a single peak for 
the two defined zones (A and B), whereas the [NII]/16584A 
emission presents a bi-modal distribution in Zone B and a 
redshifted peak in Zone A. This behaviour can indicate that the 
2+ appears to be distributed, forming a rather filled expanding 
bubble. 



Analysing the kinematics, abundances, and physical proper- 
ties, we propose an origin for this shell that could explain why it 
is only observed in the direction of the minor axis of the ellipse 
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defined above. We argue t hat ejections of m aterial in the RSG 
stage were non-spherical ( Hump hrey s|20 10 1, forming an ellip- 
soidal RSG shell denser towards its major axis. When the WR 
wind collides with the RSG shell, along the minor axis, the colli- 
sion does not find strong resistance, fragmenting it and breaking 
up the MS-surrounded bubble (as predicted in GS96 models). 
The shock then propagates in a tenuous medium where the high 
temperature favours the [OIII]/15007A emission seen and where 
a lower cooling rate could prevent the material from the WR 
stage to be observed in the optical. This idea is basically consis- 
tent with previous findings of Gruendl et al. ( 2000 1. 



5.3. Surrounding skin (SS) 

Finally, we propose the existence of an external and faint shell 
around the whole nebula, the surrounding skin. This shell is 
barely noticeable after close inspection of the INT-WFC im- 
ages, but analysis of the integrated spectra E2 and MB3 shows 
that a faint emission exists beyond the outer shell. Although the 
zones where both spectra were extracted do not present any ap- 
parent emission in Fig. [9] non negligible emission-line fluxes 
were measured on the spectra. 

Properties obtained from these two zones reveal that they 
belong to a different region than the rest of integrated spectra. 
Both appear to have similar electron temperatures and densities, 
and we find here the lowest N/H abundance (lower than solar). 
Neither helium nor neon lines were detected in these integrated 
spectra. The N/O and O/H in these regions are in the range 
of the expected values of the ISM at R G = 10 kp (NGC 6888 
galactocentric radius) considering the chemical abundance gra- 
dients. Also, I(Hy8) has its lowest intensity in E2 (except for the 
bullet, as explained below). We checked that the emission in this 
shell is intense enough not to be confused with general diffuse 
ISM emission, but it is not strong enough to belong to the OS. 
The derived oxygen and nitrogen abundances appear consis- 
tent with the ISM expected values at this galactocentric distance. 

Following GS96, we think that this skin should represent the 
early interaction between the bubble created by winds from the 
MS star with the circumstellar medium. 



5.4. The bullet 

As we explained in Sect. [T] observations in the pointing called 
Bullet were made to shed light on the nature of a dark object near 
the centre of the nebula. Measures in this zone (integrated Bl) 
indicate that it is the region where the HB flux (and other emis- 
sion lines) is lowest, but still not negligible. In this region, abun- 
dances of N and O are slightly lower than for the solar values of 
Asplund et al. ( 2009 1, but in the same range as the outer shell. 
The physical properties inferred (T e and n e ) are not remarkable 
in relation to other regions of the nebula. 

Attending to the measured intensities, we propose that 
this object is not behind the nebula because in that case all 
the nebular emission would be measured in front of it, so it is 
probably located between the main body of NGC 6888 and us, 
but not out of the nebula, because in that case all the emission 
lines would be blocked. With this analysis and considering the 
physical properties of Bl, we suggest that the bullet is placed 
in the surrounding skin defined above, and it should only be 
blocking light coming from the internal shells (see Fig. 15 i. 




Fig. 15. Sketch of the effect of the bullet on the nebular light in 
our line of sight. If the bullet is placed in the surrounding skin, it 
would block the emission from the IS and OS, but not from the 
edge of the SS. 



The observed bubble can also be explained using the GS96 
scheme. Maybe such a clump formed in the early collision be- 
tween MS winds and ISM, or perhaps as an accumulation of 
molecules and dust of the very cool RSG wind, later impelled 
by the WR shock. In any case, new observations in other wave- 
lengths, such as radio and infrared, would be necessary to under- 
stand the nature and composition of this object. 

6. Summary and conclusions 

In this article, we have presented an exhaustive study of the 
wind-blown bubble NGC 6888 with integral field spectroscopic 
observations carried out with PPAK. The most relevant results 
presented in this paper can be summarized as follows. 

1) We obtained interpolated images of four regions of the 
nebula in several wavelengths. We observe two different 
patterns for the low and high ionization ions. Some features 
are only observed in the [OIH]/15007A emission line, 
indicating clear differences in the excitation conditions. 

2) In the X-ray emitting zone of the nebula, a dithering scheme 
was adopted in the observations. This allowed us to generate 
cubes to carry out a thorough bi -dimensional analysis of 
this region: 

• We obtained a map for the reddening coefficient, c(H/3), 
finding a clear non-uniform structure, rather filamentary 
with a mean value of ~0.57. 

• We derived the electron density map, n e , showing a dis- 
tribution with a peak density (~400 cirT 3 ) spatially co- 
incident with the zone with the highest emission-line 
fluxes in the nebula. 

• We generated radial velocity maps. The velocity field de- 
rived from the Ha emission line is not spatially homoge- 
neous and is moving faster in its southwest part. 

• We studied the ionization conditions plotting pixels from 
intensity line ratios maps in three diagnostic diagrams. 
From these, we generated a map defining two main spa- 
tially well delimited zones: Zone A showing correlations 
between lines as could be expected from ionization struc- 
ture, given the different ionization degree of the ions, and 
where the maximum emission in low excitations species 
was found; and Zone B where the relations between the 
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line ratios show significant scatter in diagnostic diagrams 
that is located in the area corresponding to the peak of the 
[Oin]/15007A emission. 

We compared the statistical frequency distributions of 
the radial velocity for different lines. The [OIII]^5007A 
histogram shows a single and centred peak in both zones 
(A and B). The [NII]/i6584A presents one redshifted 
peak in Zone A and a bimodal distribution in Zone B. 
This result could be explained as the two components of 
an expanding and broken shell. 

We evaluated the possible presence of shocks to explain 
the observed behaviours in Zone B, while resorting to 
the models presented by Allen et al. (2008). We found 



that our data of Zone B can be represented well by 
models in which shock and a precursor are involved, 
assuming n = 1000 cm 4 , twice solar abundances, and 
shock velocities from 250 to 400 km s _1 . 

3) We generated nine integrated spectra extracted from the 
four pointings of the nebula. We measured line intensities to 
obtain physical parameters and chemical abundances: 

• We derived electron temperatures, T e , using the Rm = 
[NII]ii6548,6584/[NII]^5755 ratio, obtaining values 
ranging from -7700 K to -10200 K. From the 
[SH]/t/l6717/6731 line ratios, electron densities n e , were 
calculated by deriving differences of up to 350 cirT 3 be- 
tween the studied zones. 

• We inferred that nearly all the zones present a slightly un- 
derabundance of oxygen with respect to the solar abun- 
dances. We derived a nitrogen abundance similar to the 
solar value in five regions; however, in other zones N/H 
appears enhanced up to a factor of 6 or even 8. Helium 
abundance presents an enrichment in most of the inte- 
grated zones, too. 

• We compared selected emission-line ratios of the 
integrated spectra with predictions from photoionization 
models on diagnostic diagrams. The values obtained 
with these models for log (N/O) agree with the abun- 
dance derived. 

4) We provided a scenario for the evolution of the central 
WN6 star to explain the features observed in NGC 6888. 
The proposed scheme consists of a shell structure, but 
allows irregular and/or broken shells with different physical, 
kinematical, and chemical properties: 

• An inner, elliptical and broken shell (IS), with a strong 
overabundance in N/H and slightly deficient in O/H. 
We argue that the IS material corresponds to the RSG 
and WR shells that were shocked and expand in a 
low-density medium. 

• An outer and spherical shell (OS) that is weaker in all 
the observed emission lines except in [OIII]/15007A. The 
helium appears enriched here, while N/H and O/H do not 
appear enhanced or do not present any slight deficiency. 
We argue that this shell is the early MS bubble broken up 
as a consequence of the collision between RSG and WR 
shells. 

• We propose the existence of an external and faint shell 
around the whole nebula, the surrounding skin (SS). 
Measured fluxes in this regions are the lowest but not 
negligible. We think that this skin should represent the 



early interaction between the bubble created by winds 
from the MS star with the interstellar medium. The abun- 
dances inferred for this component appear to be typical 
of the local ISM. 

We checked that the dark region (bullet) observed in 
the centre of NGC 6888 appears to be associated with 
the nebula and is probably located in its most external 
shell. Observations in other wavelengths, such as radio 
and infrared, would be necessary to fully understand the 
nature of this object, but we speculate that it could be 
a dense clump impelled by the WR winds crossing the 
shells like a bullet. 
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